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Abstract 

We evaluate non factorizable 0{aas) corrections to the process 
Z — > 66 due to the virtual t-quark. All two-loop vertex diagrams with 
W^s and charged ghosts $'s are included. They are evaluated in the 
large top-mass expansion up to the 10*^ order. Gluon Bremsstrahlung 
is taken into account by integrating over the whole phase space. All 
calculations, including Bremsstrahlung, are done in dimensional regu- 
larization. The expansion coefficients of the large mass expansion are 
given in closed form. Their expansion in y = m\/Am'^ is in agree- 
ment with the coefficients up to 0{m^r /m^) as given by Har lander et 
al. i. 
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The large statistics of the LEP I experiments yielded precise values for 
the partial Z-decay width into 6-quarks; expressed in terms of the ratio 
Rb = 0.21656 ± 0.00074 § this is a precision of ~ 0.3%. Therefore precise 
high order calculations in the Standard Model (SM) are needed. 

The partial width of the Z-boson decay into a quark-antiquark pair can 
be parametrized as 

where 

f = ^("^ + <'^)(l + n',(m|))-\ (2) 

(|l]) defines the non factorizable contribution A. s, c are the sinus and cosinus 
of the weak mixing angle; iVc = 3 is the color factor; v and a are vector and 
axial couplings related to isospin I3 and charge Q of the fermion by 

v = l3- 2Qs\ a = h. (3) 

The renormalized self-energy Tlz in @) accounts for the universal correction 
due to Z-boson renormalization. 

The electroweak correction 5ew has been calculated long ago [Q] (see also 
[^]). The QCD correction factor (1 + — + . . .) is known now up to 0{al) 
in the massless limit ||^ and also at 0{mg/s) [^]. However, electroweak and 
QCD corrections do not factorize exactly. Therefore the correction factor A 
appears. This correction for the quarks of u-,d-,c- and s-type was found in 
[0], where the masses of quarks were neglected compared to m^. For the 
b-quark case it is necessary to consider diagrams with the virtual t-quark. 
The leading term proportional to has been calculated in ^ and the term 
proportional to logmf in 0. See also recent work []TU[ where leading mt 
corrections 0(m^Gi?Q;^) were calculated for Z ^ bb decay mode. 

In the following we are interested only in contributions with the top quark. 
However, this alone is neither finite nor gauge invariant (see e.g. the dis- 
cussion in |l|, 0). Instead we have to consider all diagrams with 14^-boson 
exchange. In other words we have the gauge invariant decomposition 

A = A^ + A^, (4) 

where A^'^ denotes contributions due to Z, ly-boson exchange, respectively. 
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Recently in the subleading terms for A;, were found. The expansion in 
1/mf has been made up to {1/mfY. Each coefficient of the expansion in turn 
is expanded in the small parameter y = m|/4m^ by a large mass expansion 
of subdiagrams. The calculations of [|l| were done by cutting the 3-loop two 
point function of the Z-boson. 

In a sense our approach is complementary to that of [|l|: we calculate 
directly the amplitude of the Z — > 66 process and integrate over the final state 
phase space. In addition we have to add the gluon Bremsstrahlung to form 
an IR finite quantity. We also use the standard large mass expansion (LME) 
technique |Tl| for both real and virtual gluon processes: 1/mf is considered 



as small parameter and the expansion in 1/mf has been performed up to 
(l/m^)^°. We do not, however, perform a further expansion of subdiagrams, 
therefore we obtain closed expressions for the expansion coefficients. 

However, after reexpansion we fully agree on the first few coefficients 
given in |1| . We also agree on the numerical estimates following from and 
from our work. A detailed comparison was given in [|12|. However in the 



present work the results are given in the parametrization (|1|) which seems to 
be very natural. 



In a recent paper |13] some scalar 2-loop vertex diagrams relevant for A;, 
were analysed. There it was found that the series in 1/mf has a rather bad 
behaviour (expansion coefficients grow like 4"'). Therefore one would assume 
that higher terms of the expansion are needed. However, for the complete 
physical quantity this is not the case. The 'dangerous' terms cancel as well as 
the auxiliary structures like the polylogarithms Li2, Lis which enter separate 
diagrams but not the sum of all contributions. 

In the course of our calculation we used a program written in FORM 



[Q. The input for the FORM procedures was generated by the C program 
DIANA |15[. 

The amplitude of the process Z{q) b{pi) + 6(^2) is given by 

M = e^^{q)T,{q,p,,p2), (5) 
where €^{q) is the Z-boson wave function and the amphtude reads 

6 — r 1 

= -i—b{pi)[j^v{q^) -7M75«(g^)J&(P2) (6) 
where e is the electric charge; b, b are wave functions of b-quarks. At the tree 
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level the couplings Vb and are given by 



Vb 



1 
2 



sin^ Qiv, 



(7) 



while in higher orders gets corrections to both real and imaginary parts. 

Evaluating the width from @, using the standard rules, we get in d = 
A — 26 dimensional space-time 



rrir. 



(4-d)/2 



aNcTUz T{d/2) 



{\v\'+\a\')., 



(8) 



where /2 is an arbitrary parameter with the dimension of a mass. It is related 
to the parameter /i of dimensional regularization by jl'^ = 47re~'^'^^/i^. 

At 1-loop level the diagrams contributing to this process are shown in 
Fig.l. They were evaluated in P, §]. 





W,0\ /W,<I> 



Figure 1: 1-loop diagrams with the t-quark contributing to the Zbb process 
of order 0{a). Diagrams of order 0{aas) are obtained by adding gluon lines. 

Mixed 0{aas) corrections are obtained from the diagrams of Fig.l by 
adding in all possible ways a gluon line. In this way we obtain 2-loop vertex 
diagrams as well as 1-loop Bremstrahlung diagrams. 

To obtain the renormalized expressions, for our purpose it is enough to 
renormalize the mass of the t-quark to order 0{as) and the wave function of 
the b-quark up to 0{aas)- 
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The renormalization constant of rrit in the on-shell scheme is given by 

Zm, = l~^cJ^ + l~^-\og^]. (9) 



The effect of the b-quark renormahzation is effectively performed as 

{v - a)R = Z^[v - a)o, (10) 

with the 6-quark renormahzation constant Z}^ (see appendix), while the right 
handed combination [v + a) remains unrenormalized. 

To make the width IR finite, in addition to the nonradiating process we 
have to add the soft contribution from the Bremsstrahlung process 

Z — .6 + 6 + ^. (11) 

Actually we will include the hard gluon emittion as well by integrating 
over the full gluon phase space. Of particular interest is the calculation of 
the Bremsstrahlung in terms of the LME. The kinematics of the (gluon-) 
Bremsstrahlung for the process (|11|) is given by 

g^Pi+P2+P3, v\=v\=v\ = ^- (12) 
Thus we have 3 invariants 

V\Vz- (13) 

We are interested in the integrated Bremsstrahlung. The corresponding 
width in d dimensions can be written as 

7687r3r(rf-2) 

dxdyx^~\\ - xY'^-^y^/^-\l - vY'^-^ IMb^I' , (14) 

where is the Bremsstrahlung amplitude and the invariants can be ex- 
pressed in terms of x and y as 

m| 

P1P2 = ^a;(l - y) , 




m 



P2P3 = ^(1-a;) 
2 

mi 

P1P3 = -rr^y ■ 
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It should be noted that neither (^) nor ( p!^ are finite separately after the 
renormalization. They both have IR divergences up to 0(l/£^). Therefore it 
is quite important to use general d-dimensional expressions and take the limit 



c? — >■ 4 only after adding up (|^) and (|1^) . Otherwise some finite contributions 
will be lost. 

We have checked by explicit calculation that indeed the result is finite in 
the sum of T and T^r and our result for reads 
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31 o 283 2 105443 623 7 
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(15) 



= — — + 2.3057 - 8.0270t - 28.0471t^ 



where Lt = log{mf/m'^), L^. = logc^ and y = m|/4m^. The only function 
that enters the answer is 

2 Jo ^/T^ ^ ' 

We do not give here higher coefficients in analytic form because of their 
complexity. Instead we give below 10 coefficients of the expansion 
numerically. For rrit = 175GeV, niyy = 80.33GeV and niz = 91.187GeV we 
obtain 

4.1878 
T 

- 39.5864t^ - 32.7842t^ - 8.7501t^ + 20.0429t^ 

+ 36.7551t^ + 37.3991t^ + 65.1923t^ + 307.874t^° , (17) 

where t = m^/m^ ~ 0.21. 

Separating the leading term from the higher order ones, we may write 
this as 

= 19.8764 - 1.0654 = 18.8109, (18) 

which tells us that the terms of higher order amount to only ~ 5% of the 
leading term. Obviously also the series in l/m^ converges quite rapidly. The 
reason, however, for the smallness of the correction is the alternation in sign 
of the various higher order contributions. 

For completeness and for comparison we also present the numbers for the 
non factorizable contributions of the lighter quarks and the Z-contribution of 
the b-quark. The result for these contributions are taken from 0. To slightly 
improve them, we performed a Fade approximation for the expansions in 
X = 1 for A^ and x = m|/m^ for A^ (however this gives only minor 
changes of order of few percents in formulae (|19D and (|20|) ). For A^ we have 

—0.796 for d,s,b ^ ^ 
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Taking into account (0), we have for A'^ 

-3.652 for u,c 
= { -3.745 for d,s (20) 
f 18.811 for b 

The above numbers demonstrate again that due to the heavy top quark the 
A^ contribution is larger by almost an order of magnitude than all the other 
contributions. The smallness of the subleading terms in 1/mf compared is 
surprising, however. This is a highly nontrivial result. It is, e.g., in contrast 
to the results found in Ref. |]l6l, where an 0(a^m^/m^) calculation of Ar 
was performed and it was found that this correction is of the same order as 
the leading 0{a'^mf/m^^) correction. 



A detailed comparison of our results with those of Ref. w is given in [12 



A Renormalization of the b-quark wave func- 
tion 

In this section we evaluate the virtual top quark contribution to the wave 
function renormalization constant Zb for the b-quark field in the on-shell 
scheme. We need this constant up to order 0{aas)- At 1-loop order two 
diagrams contribute. They are shown in Fig. 2. Adding in all possible ways 
one gluon line we get 8 diagrams of order 0{aas). 




b t b 



Figure 2: 1-loop diagrams with the t-quark contributing to the b-quark 
wave function renormalization in 0{a). Diagrams in 0{aas) are obtained by 
adding one gluon line. 
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The self-energy of the b-quark reads 

m = P ^l{p) +P^-Y^ ^RiP) ■ (21) 

By cxpUcit calculation we find = 0. Therefore the renormalization con- 
stant Zb in the on-shell scheme (mb = 0) reads 

Z.= l + i^S.(0) = l + i^£j(., + ^C...). (22) 

For the sake of completeness we give below the first 4 coefficients of the 
expansion of Zi^2 in m^r/m^. Note that we have to keep terms up to 0{e) in 
the 1-loop part and terms up to 0{1) in the 2- loop part. 

Zi = 

mf r 1 1 1 3 

+e(--C, - -L,K + -Lt - -a + -L, - -Ll - -)\ 
^ V 8^ A ^ " 9> 8*8'^ 8 16^ J 



1 r Ir 1 

h I/J + -L„ 



m?\4* 4 V4*^ 8 8* 4'' 8 
mtvfb 3 / 5 9 5 2,3 5 



ml 



Z2 

mf i 3 1/3. 3, 5 

-Ln + - 
2^2 



+ 2^* 

+4 + ^C2 + 3L,L^ - 5L, + ^L,2 - 5L^ + ^Lj | 



8 



, f 3 1 3 3 \ ml. 27 21 33 1 

where Lt = log(m^/m^) and = log(/iV"^vK)- 
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